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RANDOM WALK STUDY OF ELECTRON MOTION IN HELIUM IN 
CROSSED ELECTROMAGNETIC FIELDS 
by Gerald W. Englert 
Lewis Research Center 

SUMMARY 

Random walk theory, previously adapted to electron motion in the presence of an 
electric field, is extended herein to include a transverse magnetic field. Application is 
made to a weakly ionized helium gas. Interactions of the electrons with the helium 
atoms are based on integral and differential experimental cross section data. Elec- 
tronic excitations and ionizations as well as elastic collisions are included. Electron 
trajectories between collisions are based on exact solution of the equations of motion 
for a Lorentz gas. 

The restrictive effect of the magnetic field on electron motion increases the number 
of collisions per walk required for an assembly of electrons in an arbitrary initial dis- 
tribution to relax to a terminal steady state condition. This limits the maximum value 

of the ratio of magnetic field strength B to background pressure p possible to study 
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with reasonable computer time. A range of B/p from 2.25x10 to 2.25x10 weber 
per newton (30 to 300 G/torr) was investigated. This covers the range of a very weak 
to a very strong influence of B/p on electron transport coefficients. 

Time of relaxation of electron energy distribution, determined by random walks, is 
approximated by a simple expression based on energy exchange between the electrons 
and the electric field. The influence of the magnetic field is to effectively reduce the 

ratio of electric field to background pressure by^ + (wt)^ (where co and are the 
cyclotron and collision frequencies, respectively). 

Electron transport coefficients as well as energy and velocity distributions were de- 
termined and compared with existing theory. Effect of the magnetic field on mean energy 
and transport coefficients is in most cases predicted within 10 percent of the random 
walk results by use of the Hall parameter ojt. The constant mean free time and the iso- 
tropic scattering approximation often used in analytical studies, are accurate within the 
scatter of the data (estimated to less than 5 percent). 



INTRODUCTION 


study of electron motion in an electric field has been actively pursued for three- 
quarters of a century (refs. 1 and 2). It was, in fact, investigation of the cathode rays 
of the gaseous discharge which led to the discovery of the electron (ref. 3). The random 
walk concept, first formulated in 1905 (ref. 4), contributed significantly to descriptions 
of stochastic particle motions as present in such gaseous conductors. Except in the 
simplest cases (refs. 5 and 6), these descriptions lead to hard-to-solve differential 
eqviations (ref. 5). On the other hand, calculation of random walks step by step permits, 
in principle, inclusion of much physical detail with little mathematical complexity. Such 
random walks have very seldom been performed; apparently because of the burdensome 
repetitious task of calciilating the motions of statistically representative numbers of test 
particles throughout an enormous number of steps. 

Appearing first in the literature is the manual effort of Yarnold (ref. 7) to determine 
an electron energy distribution at low electric field strength. With electron-molecule 
energy exchange in mind, a low mass ratio M/m of 100 was used in an effort to reduce 
the number of operations (collisions) required in the walks. Later, with the aid of a 
computer, Wannier (ref. 8) calculated the walk of one ion for 10 000 steps, yielding the 
relatively easy to determine drift and diffusion motion. Some approximate ion velocity 
distributions were attempted. Only elastic collisions were considered. Relatively large 
samples of test particles and/or numbers of steps per walk are required to find detailed 
electron energy distributions when inelastic losses such as electronic excitation and 
ionization are considered. It is especially laborious to determine the Townsend first 
ionization coefficient at moderate ratios of electric field strength to background pres- 
sure since an electron does not have sufficient energy to cause an ionization until it dis- 
perses into the tail of the energy distribution. 

It was observed, in the recent random walk (RW) effort of reference 9, however, 
that relaxation times (thus number of steps) are inversely dependent on the coefficient of 
dispersion (second moment) of the energy exchange between the electrons and the electric 
field rather than the much smaller energy exchange between the electrons and back- 
ground molecviles. It was, in turn, demonstrated that the modern computer is capable 
of making detailed RW studies of the motion in a gaseous conductor with reasonable ex- 
penditures of computer time. Calculations for a typical set of conditions require about 
5 to 15 minutes on a modern electronic computer. 

The RW was applied to the electron motion in a weakly ionized helium gas in refer- 
ence 9 over a wide range of E/p. (Symbols are defined in appendix A.) In the present 
study this effort is extended to include a magnetic field B perpendicular to the electric 
field. The ratio B/p is varied from low values where it has a negligible effect on the 
motion to high values where the electrons are closely tied to the magnetic field lines. 
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The ratio E/p was held constant at a value of 22.5 volt-meter per newton (30 V/cm- 
torr), which is high enough for appreciable occurrences of inelastic collisions, yet well 
below the range of electron runaway. Electron drift, diffusion, and Townsend first 
ionization coefficients are determined as well as mean energies, and velocity and energy 
distributions. Pertinent theory and its simplifying approximations are appraised. 


ANALYSIS 

Each test electron of a large sample is walked for predetermined macroscopic time 
intervals. Each such walk comprises a large number of steps. Each step includes an 
analytically determined trajectory terminated by an encounter with another particle. 
Time between collisions as well as the type of interaction and concomitant scattering 
angles are determined through selection of random numbers. The macroscopic time 
intervals are selected to suitably study the relaxation of an ensemble of electrons from 
an initial arbitrary, to a steady state, velocity distribution. 

In the crossed field configuration of interest (fig. 1) both the electric field and mag- 
netic field are constant in space and time. The background helium gas is assumed only 
weakly ionized so that interactions between charged particles can be neglected. The 
electrons are quite ineffective in heating the relatively heavy background atoms. Thus 
these atoms are essentially cold, or stationary, and the medium is a Lorentz gas 
(ref. 10) . The time during a short range interaction, as between an electron and a neu- 
tral, is negligible for the gas densities of interest. Thus the effect of the electric and 
magnetic fields on the electron motion need be accounted for only during the free time t 
between collisions. 


Electron Trajectories During Free Time 

Tlie conservation of momentum of an electron during the free time is 

m — = -e(E + V X B) 
dt 

where e = 1.602X10"^® coidomb. The solution to this equation (ref. 10) for the velocity 
components is 
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= (V . + — )cos(wt) + . Sin(o)t) - — 

^ B/ B 

Vy = Vy^j cos(cot) - sin(a>t) 


v„ = v„ . 


and for the spacial coordinates is 


X 


' '‘i " w {fx. i * l) * ''y , i [' - 

*' ° ’'* * ip"* ■ (''x.i * l) 






( 1 ) 


( 2 ) 


Z = z. + V .t 
1 Z,1 


J 


where v = (v .,v .,v .) at t = 0 and, for the crossed field configuration of figure 1, 

Aji JjJ- ^,1 

E = (0,-E,0) and B (0,0,B). The subscript i denotes conditions at the start of the 
free trajectory, which are the conditions resulting from the preceding collision. Equa- 
tions (1) and (2) describe trochoidal electron trajectories having cyclotron motion of 
frequency w = -eB/m and drift velocity equal to E x B/B^. 


Free Time Between Collisions 


Descriptive integrals and approximate solutions . - Diirii^ each step of a random 
walk, a free time t is selected at random from a distribution of free times f(t) . In its 
general form this distribution can be expressed (ref. 9) as 



( 3 ) 
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where Q is the sum of the total cross sections of all the various type interactions be- 
tween the test electrons and background atoms. 

The mean free time between collisions is 


T = 




tnQ^v exp 



which reduces to (ref. 9) 



(4) 


by use of integration of parts and L'Hospital’s rule. 

The interaction energy of an electron (or ion) with a neutral atom is sometimes ap- 
proximated by the polarization potential (ref. 8) which varies as the inverse fourth 
power of the distance between centers of the two interacting particles. This is the case 
of the Maxwell molecule (refs. 10 and 11) for which Qg^v is a constant, being the 
total elastic scattering cross section. It is thus plausible that Q„v is close to being 

a 

constant, as assumed in reference 12, for example. 

A plot of the experimentally obtained value of Q v for helium (refs. 13 and 14) 

^ -14 

against e is shown in figure 2. A value of Q^v of 7. 8x10 cubic meter per second 

A 

is within 5 percent of this cxirve for 5 e ^ 80 eV; however, for e « 5 eV this is a 
poor approximation. 

When Q V is constant over the time period, equation (4) reduces to 

A 

and a distribution of free times becomes simply 

^(7) ' ““p 
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The free time random number correspondence (ref. 15) 



integrates in this case to 



The random numbers are drawn from a set of numbers uniformly distributed over the 
interval from 0 to 1. Routines to provide such numbers are available in most computer 
libraries. 

Numerical evaluation of free time integrals . - It was found in the study of refer- 
ence 9 that the components of drift velocity and diffusion coefficient parallel to the elec- 
tric field, Vj^ y and D^., were the most sensitive of the calculated electron transport 
coefficients to approximations of equations (3) and (4). Use of equation (5) gave the de- 
sired accuracy of v^ ^ and D^, only for E/p < 7. 5 volt-meter per newton (10 V/cm- 
torr). The electron energy distributions and remaining electron transport coefficients 
were found sensitive to the approximation of equation (5) for E/p values of 45 volt- 
meter per newton and greater. Results were found to be much less affected by the ap- 
proximation of equation (6). 

The presence of a magnetic field restrains electron motion and, as shown later, 
lowers the effective E/p making it less than 22. 5 volt-meter per newton (30 V/cm-torr) 
for the present study. Nevertheless the suitability of equations (5) and (6) for use in the 
present study is further appraised. Typical results showing the influence of the crossed 
field configuration on the mean free time parameter T = nQ .v.t and on the distribu- 

O. ^ 1 1 

tion of free times about the mean free time are shown in figures 3 and 4, respectively. 

The solid lines of figure 3 present T values obtained by numerically integrating 
equation (4) with use of equations (1) and the data of figure 2. 

Empirical curves, shown by the dashed lines, were fit to the numerical results for 
use in the computer program. The empirical relations used are given by equations (33) 
and (34) in table I. Figure 3(a) shows the effect of cos 8, and e. on T at (p. = 0 and 

C 111 

B/p = 7. ixlO"^ weber per newton whereas figures 10(b) and (c) show the influence of cp 
on T at various and cos 6^. Figure 10(d) shows the effect of B/p on T at var- 
ious and at the polar angle n where the effect is most pronounced. The parameter 
T was found to differ appreciably from 1 only when is less than 5 eV and B/p is 
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much less than 10 weber per newton. The electron energy distributions, to be pre- 
sented later , show that only a relatively small number of electrons lie in the low energy 
range less than 1 eV. 

A value of t/r is determined by selection of a random number for each step of 
a random walk. The corresponding value of t is then obtained by use of 


The empirical relation between t/r and R^, which approximates numerical solution of 
equations (7) and (3), and which is used in the computer program, is given in equa- 
tion (35) of table I. The distribution of t/r, obtained by drawing values of R^ from a 

uniformly distributed set and then using this empirical relation, is shown in figure 4 by 

4 

the data symbols. These symbols represent tallies of 3x10 discrete values of t/r. 

The dash-dot lines represent the distributions obtained directly by numerical integra- 
tion of equation (3) . 

The approximate distribution of equation (6) is given by the solid lines. The largest 
differences between equation (6) and the exact numerical results of equation (3) is in the 
vicinity of an value of 0.5. Equation (6) provides a good approximation to equa- 
tion (3) for ^ 5 eV. 


Type of Interaction 

At the end of each free time period a random number Rq is drawn to determine the 
type of interaction. 

It can be seen from the data of references 16 and 17 that the number of excited 
helium atoms is negligible for background gas temperatures of interest. The electron 
affinity for neutral helium atoms is less than zero, thus the formation of negative atoms 
is also negligible (ref. 18). The interactions of interest are, therefore, between free 
electrons and ground state atoms, and 

*^a = %l ^ex + *^ion 

The correspondence between the type of interaction and Rq can then be as follows. 
When 



% 
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the electron undergoes an elastic collision. When 


Q 


Q 


eZ _ D ^ 

,Rq< 


Q. 


the calculation for an excitation is made, and when 


%l '*' %x 

Q. 


R, 


Q 


an ionization is assumed. 

The integral cross sections used in the analyses are shown in figme 5. The ioni- 
zation cross section is from the experimental investigation of reference 19. The exci- 
tation cross section is the sum of the cross sections for excitation from the ground state 
to the next seven energy levels in both the singlet and triplet spectral series of helium 
(ref. 17). The cross section for elastic collisions was obtained by subtractii^ Q 
and from Q^. 


Scattering Angles 

Once the type of interaction is known, the scattering angle can be determined. The 

relation between differential scattering cross section a(e, x) and uniformly distributed 

random number R is 
X 




/•X 

/ a(e, x')sin X’ dx’ 



rif 

I a(e,x')sin x' dx’ 

n 


where x is the scatter ir^ ar^le in local spherical coordinates. 

The experimental differential scattering cross section data of references 20, 13, 
and 21 were used for elastic encounters (fig. 6(a)) and that of references 20, 22, and 23 
for inelastic encounters (fig. 6(b)). 

The inelastic data are primarily the cross section for excitation of the most prob- 
able state (2^P) from the ground state (l^S). Little differential scattering data to other 
states is available. The data of reference 24 show essentially the same distribution of 
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scattering angle for excitation to the 2^P, 3^P, or 4^P states. The data of figure 6(b) 
are also used to represent ionization events. 

The recent experimental restilts of reference 23 show a steeper rise of ct(x) with 
decrease of x and a slightly lesser effect of e than the results of references 20 and 22. 
Two sets of curves were therefore used to find the effect of these differences on electron 

4 

transport coefficients. The data symbols in figure 6 represent tallies of 3X10 values of 
X obtained from drawing R and using relations (36) to (38) of table 1. 

A. 

Azimuthal angle a is uniformly distributed; thus 


Q! = 2T7R 

a 


( 10 ) 


Energy Losses 

The energy lost from an electron Ae and given to an atom during an elastic inter- 
action is (ref. 10) 


Ae = ^ e(l - cos x) (11) 

M 

For excitations the average Ae for scattering to the first seven states in the singlet and 
triplet series is used. This average was fo\ind equal to 21.9±0. 1 eV over the range of 
interest (ref. 9). 

The minimum energy lost in an ionization event is the potential energy, 24.46 eV. 

It is assumed that the remaining kinetic energy is evenly divided between the test and 
the liberated electron. Thus, = (e - 24.46)/2 after an ionization event. The effect 
of various ionization energy loss assumptions on RW results are studied in reference 9. 


Coordinate Systems 

The transformation to spherical coordinates at the end of each free time period is 
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(12) 


-iAx\ 

cp = tan ( — 
Vz/ 


whereas the transformation from the scatterir^ coordinates back to spherical laboratory 
coordinates at the end of each collision calculation is 






v2 -lAl 


m 


0. = cos“^ ^cos X cos 6 + sin x sin 6 cos aj 




<p^ = (p + 


cos"^ ^COS X - cos 9 cos 0.^ 


sin 0 sin 0. 




(13) 


In summary, to perform a random walk for this weakly ionized plasma, four ran- 
dom numbers are drawn for each step. The first, R^, to determine the ratio of free to 
mean free time, then Rq to determine the type of encounter; finally R^ and R^ to 
determine the scattering and azimuthal angles. The trajectories durii^ the free time 
are determined by analytical solution to the equations of motion, whereas the collisional 
interactions are based on experimental data. 

As an initial condition at the start of the very first free path for each test electron, 
an isotropic distribution in angle 0^ was assumed. The electron energy at the start of 
each walk was set equal to a constant value close to the expected average energy at the 
steady state. The steady state, or terminal conditions were found in preliminary cal- 
culations to be independent of the initial condition selected. ^ 
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Determination of Macroscopic Quantities at Steady State Conditions 


During the walk of each test electron, the time T equal to the accumulative free 
times of the walk, is monitored. At prescribed time periods, T = Tg, . • . , T., 
the energy distribution of all the test electrons is determined. This is obtained by 
tallyir^ the numbers of test electrons which have energy within the various equal in- 
crements of the rai^e of e. The time period was increased logarithmically until there 
were no observable changes in the energy distribution F^(e) with time, indicating that 
steady state conditions are reached. 

After such a terminal state is reached, the drift velocities can be obtained from 


^D.x 


_ - =^1' _ E 






A(T, 




and y (14) 

_ Z 

""^’y A(T. - T._j) 

where the subscript j denotes the time period, A is the number of test electrons, 
and the summations are over the free times of the A electrons between the j and j-1 
time periods. 

In like manner the various components of the diffusion coefficient are obtained from 
the second moments of distance during the free time periods as 


D 1 

^ 6 A(T. - T._j) 

y 6 A(T. - T._^) 

o X E(--i)^ 

^ 6 A(T. - T._p 


(15) 
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The number of ionizations are also recorded during the time periods so that the 
Townsend ionization coefficient Q!,p can be obtained from 


number ionization events 


(16) 


Mean energy e is obtained from 


e = 



(17) 


The walks of at least 1000 test electrons were determined for each B/p and physi- 
cal model of particle behavior considered, from which the macroscopic quantities of 
equations (14) to (17) were obtained. 


RESULTS AND DISCUSSION 
Electron Energy and Velocity Distributions 


The relaxation of the electron energy distributions from delta functions to terminal 
steady state conditions is shown in figure 7 for a wide range of B/p values. The time 
parameter tp and average number of collisions experienced per test particle in re- 
laxing to the various intermediate distributions are also listed. 

The values of these variables for relaxation to terminal conditions compare well 
with those obtained by use of the following theoretical relations based on the energy ex- 
change criterion of reference 9; 


N 


e 




(18) 


and 
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where n/p is the conventional ratio of number density to pressure at 0^ C, equal to 
2.65x10^® newton"^ meter"^ (3. 54x10^® cm“^ torr"^) and cot is the Hall parameter. 
This theory is based on the second moment of the energy exchange between the electrons 
and the electric field. The extension to include a magnetic field for the case of interest 
is given in appendix B. Note that the magnetic field serves only to reduce the E/p 

ratio by + (cot)^. The product p'^1 + (cot)^ may be considered an effective pres- 
sure (ref. 25). 

In relaxing to the steady state, the peaked distributions spread out very rapidly at 

first. At high B/p there is a shifting of the energy distribution toward lower energies 

near the final stages of the relaxation. This is apparently due to the decreased net 

amount of energy the electrons receive from the electric field as B is increased (see 

eqs. (B3) and (B4) of appendix B) which, at steady state conditions, must balance the 
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electron energy lost in collisions. At a B/p of 2x10 weber per newton, the number 
of collisions for relaxation to a terminal distribution is about five times that required 
at B = 0. Beyond a B/p of 2X10"^ weber per newton, the required number of steps 
per walk and thus computer time rapidly increased to excessive values. 

Steady state energy distributions were fovind to be independent of the energy of 
the initial delta function distributions. That the average number of steps per walk re- 
quired for relaxation was also quite independent of is illustrated in figure 8. 

Marginal distributions of the three rectangular components of random velocity, de- 
fined as 


CO - oo 


1 


f(v V V )dv dv 
' x,r y,r z,r^ y,r z,r 


oO -oo 


. / / 


^(^x,r^,r^z,P'^^x,r ‘i^z,r 




( 20 ) 


and 


oO /• 


f(v V V )dv dv 
^^''x,r y,r z,r^ y,r ''z,r 
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are plotted in figure 9 at a high and a low B/p. The distributions are close to 
Maxwellian (referenced to = 3 !<0/2) at all B/p investigated. This agrees with the 
theory of reference 10 for Qg^v equal to a constant. 


Mean Random Energy and Transport Coefficients in Helium 


The mean random electron energy and Townsend's first ionization coefficient at 
steady state conditions are plotted against B/p in figure 10. The dashed lines were 
obtained by reading el and ar^/p from the B = 0 curves of reference 9 at effective 


ratios of electric field to pressure equal to E 


/l-f 


The effect of B/p on 


~e and ar^/p is quite well accounted for by such a procedure, which is consistent with 
the findings of reference 25 for hydrogen. The mean energy predicted by reference 10, 
however, is about a factor of 10 too high. 

Two simplifying assumptions, often used theoretically, are also appraised; the use 
of an isotropic distribution of scattering angle x and the assumption that Q„v is a con- 

a 

stant (constant mean free time between collisions). These assumptions changed the 
mean energy less than 10 percent; tending to lower at low B/p and raise at 
high B/p. The Or^/p results were also charged less than 10 percent except at the 
highest B/p value. 

The oirp/p values appear to have the most data scatter of the RW results. This is 
due to the relatively small number of ionizing encounters. 

Drift velocities in the E and ExB directions are plotted in figure 11. The effect 
of B/p is shown to be quite well predicted by the theoretical relations (ref. 10) 


^D,x = ^D,x(B = 0) 


COT 




+ (cot)‘ 


( 21 ) 


and 


V 


D,y 


^1 + (cor)^ 


( 22 ) 


which are based on the assumption that Q v is a constant. Drift velocity at B = 0 was 

3 . 

obtained from reference 9. 

The use of the simplifying assumption that Q v is constant had very little effect on 
the RW determined drift, the changes being less than 10 percent. This assumption 
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raised Vp, slightly. The isotropic scattering approximation also consistently low- 
ered Vp ^ and raised This is most pronounced at low B/p and is due to the 

fact that the isotropic approximation destroys the persistence of velocity trend exhibited 
in the differential scattering cross section data (fig. 5). 

The three cartesian components nf diffusion coefficient are plotted in figure 12. 
Theory based on Q„v equal to a constant yields (ref. 26) 

a 


P°x = P^x(® = — (23) 

1 + (wt)^ 

pD (B = 0) 

pDy = (24) 

1 + (cot)^ 

and 

pD^ = pD^(B = 0) (25) 

Values of diffusion coefficient at B = 0 were obtained from reference 9. These rela- 
tions slightly underestimate the decrease of and with increase of B/p as 
shoAvn by the line comprised of long dashes. This discrepancy is more apparent for the 
z component and is due to the decrease of T with B/p (see fig. 9(a)) as yet unaccounted 
for. 

2 

Diffusion coefficient is proportional to v ; thus multiplying the theoretical values 
in equations (24) and (25) by ~ 0) reduced the analytical results to close agree- 

ment with the basic RW calculations for all three components. This is shown by the line 
comprised of alternating long and short dashes in figure 11. 

Results were very insensitive to the isotropic scattering assumption as well as to 
the use of Q v set equal to a constant. Results agree to within the amount of scatter 
in the data, estimated to be less than 5 percent. 

In general, holding Q„ constant (equal to 7. 8xl0~^^ m^/sec) over the whole range 
of € to which the electrons have access has a fairly small influence on transport coef- 
ficients. The simplification of holding Q^v constant (equal to Q„ -v.) only during any 
given free time period, and thus permitting the use of equations (5) and (6), should be of 
even smaller consequence. The effect on the transport coefficients of using equation (5) 
in place of equation (4) is observable at lower B/p, however very little effect is ob- 
servable at the intermediate values of B/p (TxlO”^ Wb/N). The change of results due 
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to replacing equation (3) with equation (6) is within the data scatter even at the lowest 
B/p studied. 

In much the same manner, the incorporation of the relatively small difference of re- 
sults between the differential scattering cross sections of reference 23 with references 
20 and 22 was found to have a negligible effect compared to the replacement of the cross 
sections of figure 6 with an isotropic distribution of scattering angle. 

CONCLUSIONS 

The present random walk study made over a wide range of magnetic field strengths 
leads to the following conclusions regarding electron motion in helium at an E/p of 
22. 5 volt-meter per newton (30 V/cm-torr) where E is the electric field strei^th and 
p is pressure at 0° C; 

1. The restraining effect that the magnetic field has on the electron motion causes 
an increase in the relaxation time and corresponding number of steps (free time periods) 
for an assembly of electrons to reach a terminal steady state distribution. This trend 
can be predicted by use of an effective ratio of electric field to background pressure, 

equal to e/ + (cot)^ (where co and are the cyclotron and collision frequency, 

respectively), in the relaxation equation for magnetic field strength B = 0. The amount 

of computer time to perform random walks sufficient to simulate the electron motion of 
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interest is within reason for B/p values up to at least 2.25x10 weber per newton 
(300 G/torr) . 

2. Distribution of random velocity was found to be close to the Maxwellian form 
theoretically predicted by Chapman and Cowling for the case of the collision frequency 
parameter Q v equal to a constant. This theory, however, predicts a mean random 
energy, about which the particles are distributed, which is a factor of about 10 too high at 
E/p = 22. 5 volt-meter per newton (30 V/cm-torr). 

3. The use of the effective pressure concept enables determination of the dependence 
of mean random energy and Townsend ionization coefficient on B/p from knowing only 
their dependence on E/p at B = 0. 

4. Existing theory (see Chapman and Cowling) was quite satisfactory in predicting 
drift motion as well as the x and y components of diffusion. The drop off of the z 
component of diffusion coefficient pD with increase of B/p could be accounted for 
through the fall off of mean random energy with increase of B/p. 

5. The simplifying assumption of constant mean free time was appraised by incor- 
porating a constant value of absorption cross section times electron velocity 

Q V = 7. 8x10"'^'^ cubic meter per second in the random walk (RW). Results in general 
were quite insensitive to this restriction and changed less than 10 percent except for the 
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value of Q!,p/p (where is the Townsend ionization coefficient) at the highest B/p 
(equal to 2.25x10”'^ Wb/N). The assumption of an isotropic distribution of scattering 
angle was studied in like manner. The largest influence of this simplification was in the 
lowing of the y component of drift velocity due to the elimination of the persistance of 
velocity in the negative E direction during scattering, 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 11, 1971, 

112 - 02 . 
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APPENDIX A 


SYMBOLS 

[The International System of Units (SI) can be used throi^hout the equations. Nu- 
merical results are often presented in other imits as well for ease of comparison with 
previous investigations. ] 

A number of test electrons or walks 
a constant defined in table I 

B magnetic field strength 

C constant defined in table I 

D diffusion coefficient 

E electric field strength 

e absolute value of electrostatic charge of an electron, 1.602x10 C 
F distribution function 

f distribution function 

-27 

M mass of a helium molecule, 6,695X10 kg 

-SI 

m mass of an electron, 9. 108x10 kg 
N average number of steps per walk 
n number density at 0° C 

p pressime at 0° C 

Q integral cross section 

R random number 

r variable defined in table I 

S variable defined in table I 

T time at which data are recorded 

t time 

V velocity 

v^g£ reference random velocity, ^2 e^/m 
X Cartesian coordinate antiparallel to E x B 
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y Cartesian coordinate antiparallel to E 

z Cartesian coordinate parallel to B 

Of azimuthal scattering angle (see fig. 1) 

Q-rj, Towsend ionization coefficient 

j3 variable in table I 

e electron energy 

© temperature 

9 polar angle 

a differential cross section 

T mean free time between collisions 

T mean free time parameter equal to ratio of mean free time to mean free time 
for constant 

<p azimuthal angle inx,y plane 

X scattering angle 

CO electron cyclotron frequency 

Subscripts: 

a absorption 

D drift 

el elastic 

ex electron excitation 

i at the start of the free time period 

ion ionization 

o condition at start of walk 

Q denotes a random number drawn to determine type of collisional interaction 

r random 

t time 

X X component 

y y component 

z z component 

a denotes a random number drawn to determine azimuthal angle 
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e 


relaxation to terminal distribution based on energy exchange 
pertains to azimuthal angle in x,y plane 

denotes a random number drawn to determine scattering angle 


<P 

X 
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APPENDIX B 


THEORETICAL APPROXIAAATION OF RELAXATION TIME 

The energy relaxation time of reference 27 is defined as 

t^=^^4r (Bl) 

(Ae)2 

In like manner the average number of collisions per test electron for relaxation to a 
steady state is 


N 


€ 



(B2) 


It was found in reference 9 that this criterion gives good comparison with RW re- 
sults if the energy exchange Ae is based on the interaction of the electrons and the 
electric field. 

Considerable energy is usually exchanged between the electron and the field during a 
free time period. That is, the absolute value of Ae is usually large. When accounting 
for algebraic sign, however, the much smaller average net amount of energy given to 
the electrons by the field just balances the energy transfer from the electrons to the cold 
neutrals at steady state conditions. An average based on the second moment of Ae 
yields the square of the gross electron field energy exchange. Relative to the magnitude 
of this average, the drift energy and energy loss to neutrals (first moments) can be ne- 
glected within the accuracy desired for relaxation time estimates. 

The effort of reference 9 is extended herein to include a magnetic field. For the 
cross field configuration of interest 


where 


(Ae)2 = (Ay)^ 


Ay = 1 


cv 


sin(cut) - 



(B3) 
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Since Ay is a function of v^. and t, 


(iy)^ 





X 




where within desired accuracy 

iM = i 

T 

as previously discussed. 

Due to the large ansrle scattering in a and x, the v . and v • are essentially 

j 1 y 5 1 

uncorrelated, thus 


f (v .V — f (v -)f (v •) 
v'' x,i y,i^ v,x'' x,i' v,y'' y,i' 


If drift velocities are neglected, f and f are even functions of v . and v ., 

VjA V,y y 7 ^ 

2 

respectively. Since the drift motion E x B/B can also be neglected for purposes at 
hand, the first moments and, therefore, also products of first moments of Ay are 
zero. Thus 


(Ay) 2 = 


n oO 

v2 . 


sin2(cot)f^^y(Vy^.)e“^/'^ dVy^. dt 


(jJ T''0 


n oo 

v2 [l - 2 cos(o;t) + cos2(cot)]f^ (v^ )e‘^/^dv^ dt 


This integrates to 


(Ay) 


O 2 2 

2 2v.t 


i + (wt)^ 


(B4) 


whei-e it is assumed that 
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V.. 


= .=iv2 


X,I y,i 


for the near isotropic distribution of interest. 
It can be fvirther approximated that 






h2q2 


(B5) 


as on page 19 of reference 26, where I is the mean free path. Thus, using equations 
(B3), (B4), and (B5) in equation (B2) yields 


N 


e 



1 + (cot)^ 


(18) 


Usir^ equations (5), (18), and (B2) in equation (Bl) gives 



-13 


Let, for helium, Q„v. = 0.78xl0~ cubic meter per second; then 

3 . 1 


10 


_ / 

N-sec/m^ 

E ■ 


(19) 
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Figure 2. - Collision frequency parameter 
vQg in helium plotted against electron 
energy e. 
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(a) Influence of electron energy and polar 
angle 9j when azimuthal angle tPi'Dand 
ratio of magnetic field strength to back- 
ground pressure B/p = 7. 1x10"^ weber per 
newton (95 G/torr). 

Figure 3. - Effect of conditions at start of 
electron trajectories on mean free time 
parameter t. 



Increment of mean free time parameter, Increment of mean free 
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Electron energy, 



(b) Influence of fj and «j when cos 0j = 0 and ratio of mag- 
netic field strength to background pressure B/p = 7. lxlO‘5 
weber per newton (95 G/torr). 



30 40 50 60 80 100 200 300 

Ratio of magnetic field strength to background pressure, 
B/p, G/torr 


(c) Influence of “Pj and cosBj when €j = 0. 5 eV and ratio of Id) Influence of B/p and when azimuthal angle Vj = nl2 

magnetic field strength to background pressure B/p= radians and cos0j = -l. 

7. 1x10'^ weber per newton (95 G/torr). 


Figure 3. - Concluded. 



Orstriiiutlon function, 





(b) Effect of 9j at electron energy «j = 
1.0 eV and cos 0; = 0. 




Free time/mean free time, t/t 

(c)Effectof 9j at «t = 0.5 eV and azimuthal (d)Effectof 9j at £j = 5.0 eV and azimuthal 

angle 9, = tt/2 radians. angle 9; = tt/2 radians. 

Figured. - Effect of conditions at start of electron trajectories on distribution F* of free time t about 
mean free time t. Ratio of magnetic field strength to background pressure Bip = 7.1xl0'^weber 
per newton (95 G/torr). Spbols denote tallies of distributions of random numbers obtained by using 
equation (35). Dash-dot lines denote numerical integration of equation (3). 
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Figure 5. - Cross sections for interaction of 
electrons with helium. Elastic cross sec- 
tions obtained from absorption cross sec- 
tions of references 13 and 14 minus ex- 
citation and ionization cross sections. 
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Scattering angle, x. deg 
(a) Elastic scattering. 





(b) Inelastic scattering from ground to 2'P state. 

Figured. - Differential scattering cross sections a versus scattering angle x- Lines are for experimental data and symbols are for tallies of 
distributions of random numbers Rx obtained with use of equations (36) to (38). 
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(a) Magnetic field strength B = 0; 
average number of collision pre- 
dicted by equation (18) Ne = 93; 
corresponding time predicted by 
equation (19) t^p = 0. 45x10'^. 


(b) Ratio of magnetic field strength 
to background pressure B/p = 

2. 25x10'^ w/eber per newton 
(30 G/torr); average number of 
collisions predicted by equation (18) 
Ne = 97; corresponding time pre- 
dicted by equation (19) t^p = 

0. 47x10'5. 


Figure 7. - Electron energy distribution at various times of relaxation from an 
initial delta function. Time t is multiplied by pressure p of helium refer- 
enced to 0° C. E/p = 22.5 volt-meter per newton (30 V/cm-torr). 



00 




Random electron energy, €[•, eV 


(c) Ratio of magnetic field strength to 
background pressure Blp = 0. 71x10“^ 
weber per newton (95 G/torr); average 
number of collisions predicted by 
equation (18) = 117; corresponding 

time predicted by equation (19) t^p = 

0. 56x10'5. 


(d) Ratio of magnetic field strength to 
background pressure B/p = 2.25x10“^ 
weber per newton (300 G/torr); average 
number of collisions predicted by 
equation (18) N« = 318; corresponding 
time predicted by equation (19) tep = 

0. 15x10'^. 


Initial electron 
energy, 



Figures. -Influence of initial electron energy £g on mean electron 
energy « versus average number of collisions per walk N. The 
theoretical value Ns for relaxation to terminal conditions is from 
equation (18). Ratio of magnetic field strength to background pres- 
sure B/p = 2. 25x10'^ weber per newton (30 G/torr). 


Figure 7. - Concluded. 
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(b) y Component of random velocity. 



Random velocity ratio, Vy |-/Vpgf 
(c)z Component of random velocity. 
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Figure 10. - Effect of magnetic field on mean random electron energy and on 
Townsend first ionization coefficient divided by helium pressure referenced 
to 0° C. 


Figure 9. - Electron marginal velocity distri- 
butions in helium. E/p= 22. 5 volt-meter 
per newton (30 V/cm-torr). 
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o Basic random walk (RW) calculation 
□ RW with iscftropic scattering 

o RW with constant QgV 

3x10^ Theory of ref. 10 
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0 


o. 


Solid symbols denote asymptotic values at B = 0 
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Figure 11. - Effect of magnetic field on drift velocities. 
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Diffusion coefficient, Diffusion coefficient, Diffusion coefficient, 



d 

Q. 
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Figure 12. -Effect of magnetic field on diffusion coefficients multiplied by 
pressure of helium referenced to 0° C. 
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